chan = sort([7 9 15 19 25 20 10 26 16 8 14 24 6]);
timerange = 501:747;
data1 = [];
data2 = [];

for index = 1:length(chan)
    data1 = [data1 mean(STUDY.changrp(chan(index)).erpdata{1,2}(timerange,:)) ];
    data2 = [data2 mean(STUDY.changrp(chan(index)).erpdata{2,2}(timerange,:)) ];
end;

[t df p] = statcond({ data1 data2 }, 'mode', 'bootstrap', 'naccu', 2000)
figure; hist2(data1,data2,20);
